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TECHNICAL  NOTE  28l4 


THE  APPLICATION  OF  PLANING  CHARACTERISTICS  TO  THE 
CALCULATION  OF  THE  WATER-LANDING  LOADS  AND 
MOTIONS  OF  SEAPLANES  OF  ARBITRARY 
CONSTANT  CROSS  SECTION 
By  Robert  F.  Smiley 


SIMMARY 


The  general  equations  governing  the  fixed-trim  water  landing  of  a 
straight-keel  seaplane  with  a  hull  of  arbitrary  constant  cross  section 
are  presented  in  such  a  form  that  the  landing  motions  and  loads  are 
expressed  in  terms  of  the  steady-planing  characteristics  of  the  sea¬ 
plane.  In  order  to  verify  the  general  validity  of  these  equations, 
solutions  are  made  for  the  water  landing  of  a  rectangular  flat  plate 
and  are  compared  with  experimental  impact  data.  Calculated  and  experi¬ 
mental  time  histories  of  draft,  velocity,  and  load  are  in  good  agreement. 
A  survey  is  made  of  the  available  information  on  seaplane  planing  char¬ 
acteristics  which  is  suitable  for  use  with  the  analysis  of  the  paper. 

INTRODUCTION 


The  National  Advisory  Committee  for  Aeronautics  has  undertaken  an 
extensive  program  of  theoretical  and  experimental  research  on  hydro- 
dynamic  impact  loads  in  order  to  establish  a  more  rational  foundation 
for  water-loading  requirements  for  the  design  of  seaplanes.  Most  of  the 
results  of  this  program  to  date  are  contained  in  references  1  to  6.  The 
development  of  the  theory  in  these  various  papers  usually  proceeds  sub¬ 
stantially  as  follows:  First,  a  theoretical  or  semiempirical  analysis 
is  made  for  the  hydrodynamic  forces  acting  during  the  two-dimensional 
impact  of  a  body  on  a  smooth  water  surface;  the  three-dimensional  impact 
and  planing  case  is  then  treated  by  assuming  that  the  fluid  flow  occurs 
primarily  in  two-dimensional  planes  oriented  normal  to  the  keel  and 
applying  an  approximate  over-all  correction  to  account  for  the  differ¬ 
ence  between  the  two-dimensional  and  three-dimensional  cases.  This  type 
of  approach  to  the  impact  and  planing  problems  has  been  found  to  provide 
fairly  reasonable  estimates  of  the  impact  loads  on  certain  types  of  sea¬ 
plane  hulls  (refs.  1  to  6),  particularly  those  with  scalloped  bottoms 
and  V-bottoms  when  the  chines  are  not  immersed  below  the  water  surface. 
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However,  for  impacts  involving  chine  immersion,  including  impacts  of  a 
rectangular  flat  plate,  accurate  two-dimensional  solutions  have  not  yet 
been  derived  for  all  cases.  For  such  cases,  or  for  cases  where  greater 
accuracy  is  desired  than  can  be  obtained  from  the  two-dimensional  analogy, 
other  procedures  for  computing  impact  loads  and  motions  must  be  developed. 
The  purpose  of  this  paper  is  to  develop  such  a  procedure  by  relating  the 
basic  seaplane  impact  equations  to  the  planing  characteristics  of  a  sea¬ 
plane  and  to  present  the  solution  of  these  equations  in  such  a  form  that 
the  impact  loads  and  motions  may  be  calculated  from  these  planing  char¬ 
acteristics.  This  derivation  begins  with  the  same  assumptions  as  were 
made  in  the  preceding  papers  on  impact  theory  (refs.  1  to  6),  namely, 
that  the  instantaneous  forces  during  a  landing  and  in  planing  depend 
only  on  the  components  of  motion  normal  to  the  keel.  The  differential 
equation  of  motion  based  on  these  assumptions  is  presented  and  then  inte¬ 
grated  to  obtain  equations  for  the  time  histories  of  the  draft,  velocity, 
and  hydrodynamic  load  during  an  impact.  The  validity  of  the  resulting 
equations  is  then  tested  by  comparisons  of  experimental  and  calculated 
impact  loads  and  motions  for  water  landings  of  a  rectangular  flat  plate. 
Finally,  a  survey,  and  evaluation  is  made  of  the  available  information 
on  the  planing  characteristics  of  seaplanes. 


SYMBOLS 


2 

(Wetted  length  at  keel) 

A  hydrodynamic  aspect  ratio,  - = - 

Wetted  area  projected  normal  to  keel 

b  beam  of  model,  ft 

Cg  function  of  angle  of  dead  rise  (0.750  for  3  =  30°) 

F^  hydrodynamic  force  normal  to  keel  (normal  to  surface  for  a  flat 

plate ) ,  lb 

Fz  vertical  hydrodynamic  force,  lb 

g  acceleration  due  to  gravity,  32.2  ft/sec^ 

H  perpendicular  distance  between  keel  and  plane  of  chines,  ft 

m  mass  of  model,  slugs 

iru  virtual  mass,  slugs 


impact  load  factor,  Fz/mg 
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V 
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(3 
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£ 

£ 

£ 

Xch 


wave-rise  factor  (see  fig.  1(c)) 

time  after  water  contact,  sec 

instantaneous  resultant  velocity  of  model,  fps 

instantaneous  velocity  of  model  parallel  to  undisturbed  water 
surface,  fps 

instantaneous  draft  of  model  normal  to  undisturbed  water 
surface,  ft 

instantaneous  velocity  of  model  normal  to  undisturbed  water 
surface,  fps 

instantaneous  acceleration  of  model  normal  to  undisturbed  water 
surface,  ft/sec^ 

angle  of  dead  rise,  deg 

instantaneous  flight-path  angle  relative  to  undisturbed  water 

surface,  tan-1  deg 
x 


perpendicular  distance  between  keel  and  point  of  initial  water 
contact,  ft 

instantaneous  velocity  of  model  normal  to  keel  (normal  to  model 
surface  for  a  flat  plate),  x  sin  t  +  z  cos  t ,  fps 

instantaneous  acceleration  of  model  normal  to  keel  (normal  to 
model  surface  for  a  flat  plate),  ft /sec2 

wetted  length  at  chines,  beams 


length  of  model  below  undisturbed  water  surface, 
beams 


z 

b  sin  T  * 


Xe  wetted  length  at  keel  or  longitudinal  center  line  of  model, 
beams 


V 


mean  wetted  length, 


^e  +  ^ch 


beams 


instantaneous  velocity  of  model  parallel  to  its  longitudinal 
center  line,  x  cos  t  -  z  sin  t,  fps 
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Dots  are  used  to  indicate  derivatives  with  respect  to  time;  for 
example, 


z  =  z  = 

dt  dt 

ANALYSIS 

Derivation  of  Equation  of  Motion 


In  the  first  part  of  the  analysis,  the  basic  differential  equation 
for  the  hydrodynamic  loads  and  motiofts  occurring  during  the  oblique 
water  impact  of  a  straight-keel  seaplane  with  a  hull  of  arbitrary  con¬ 
stant  cross  section  is  presented  and  then  converted  into  such  a  form 
that  the  equation  is  expressed  in  terms  of  the  conventional  planing 
coefficient  C-g.  Throughout  this  analysis  the  seaplane  is  assumed  to 

remain  at  a  fixed  trim  and  to  have  zero  roll  and  yaw.  The  wing  lift 
force  acting  on  the  seaplane  is  assumed  to  be  equal  to  the  weight  of 
the  seaplane,  so  that  the  net  force  acting  to  accelerate  the  seaplane 
is  the  hydrodynamic  force. 

During  the  landing  of  a  long  narrow  body  of  constant  cross  section 
on  a  smooth  water  surface,  the  hydrodynamic  force  on  the  body  is 

generally  assumed  to  be  composed  of  two  terms,  one  proportional  to. the 
square  of  the  component  of  the  model  velocity  normal  to  the  keel  £  (see 
fig.  1)  and  one  proportional  to  the  normal  deceleration  of  the  body  £ 
(see,  for  example,  ref.  l),  or 

Ff.  =  bI2  +  E£  (1) 


where  D  and  E  are  coefficients  that  depend  on  the  instantaneous 
gometrical  conditions  (body  shape,  trim,  and  draft).  The  last  term  of 
equation  (l)  can  be  considered  to  represent  the  inertia  force  corre¬ 
sponding  to  the  acceleration  of  the  virtual  mass  of  water  associated 
with  the  impacting  body;  the  coefficient  E  is  thus  equal  to  the  virtual 
mass  of  water  associated  with  the  body  (E  =  mw).  A  procedure  for  com¬ 
puting  this  quantity  is  given  in  appendix  A.  The  coefficient  D  can  be 
interpreted  in  terms  of  the  conventional  planing  coefficient  Cg  as 

follows:  In  steady-planing  tests  the  deceleration  £  is  zero  and  the 

normal  velocity  is  related  to  the  planing  velocity  Vpl  by  the  rela¬ 
tion  £  =  Vpl  sin  t  (see  fig.  1(a)).  Also,  the  normal  hydrodynamic 
force  (neglecting  friction)  is  related  to  the  vertical  force  by  the 
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relation  Fz  =  Fj.  cos  t. 
gives  for  planing 


Substituting  these  relations  into  equation  (l) 


P  P 

Fz  =  DVp^  sin  t  cos  t 


Since  the  planing  coefficient  Cg  is  defined  as 


p  _  _ £ _ 

B  \  pvPi2b2 


(2) 


(3) 


a  comparison  of  equations  (2)  and  (3)  shows  that  D  is  related  to  Cg 
by  the  equation 


D 


,2 

pb  Cg 


2  sin^T  cos  t 


(4) 


By  use  of  this  relation  (and  E  =  mw),  equation  (l)  can  be  rewritten  as 

2  2 

pb  , 

F  =  - (5) 

b  2  sin  t  cos  t 

Since  the  wing  lift  on  the  model  has  been  assumed  to  be  equal  to  the 
weight  of  the  model,  the  hydrodynamic  force  is  the  net  force  acting  on 
the  model,  and  according  to  Newton's  second  law  this  force  must  be  equal 
to  the  inertia  reaction,  or 

F  *  =  -ml  =  - -  1  +  (6) 

^  2  sin  t  cos  t 


Equation  (6)  is  the  general  differential  equation  for  the  force  and 
motions  for  a  seaplane  landing,  expressed  in  terms  of  the  conventional 
planing  coefficient  Cg. 
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Solution  of  Equation  of  Motion 

In  order  to  solve  the  differential  equation  of  motion  it  is  con¬ 
venient  first  to  rewrite  equation  (6)  in  the  form 


F  f  =  -m£  = 


2  s  in^T  cos  t  ^1  + 


Then,  converting  from  normal  force  to  vertical  force  by  multiplying 
through  by  cos  t  and  abbreviating 


“w 

1  +  — 
m 


gives  the  relation 


Pi  ' 

Pb  C* 

-m£  cos  t  =  - £ 

2  sin  t 


The  quantity  Cg  for  a  given  hull  shape,  weight,  and  trim  depends 
on  only  one  variable,  the  draft.  On  the  other  hand,  £  and  £  depend 
on  the  two  variables  £  and  tj  equation  (8)  is  thus  an  equation  in 
three  unknowns.  In  order  to  solve  this  equation  one  of  these  variables 
must  be  expressed  in  terms  of  the  other  two.  The  variable  which  is 
simplest  to  eliminate  is  £,  which  can  be  eliminated  by  the  following 
substitutions  (see  fig.  1(b)): 

• 

£  =  — - —  +  |  tan  t  (q) 

cos  T 


where  |  is  the  component  of  the  seaplane  velocity  parallel  to  the  keel, 
which,  from  the  previous  assumptions  of  no  friction  drag  and  equal  wing 
lift  and  float  weight,  is  assumed  to  remain  constant  throughout  any 
landing.  Equation  (8)  can  then  be  expressed  in  the  form 


+  I  tan  t 


2  i 

pb^CB 

2~ 

2m  sin  t 
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pb2CB' 


- . - —  u.^  -  2  2 

(z  +  |  sin  t)  2m  sin  t  cos  t- 


Since  z  dz  = 
in  the  form 


—  dz  =  —  dz  =  z  dz.  equation  (ll)  can  he  also  express 

at  at 


z  dz  =  Pb2cBtdz 

(z  +  i  sin  t)  2ms  in  t  cos  t 


-2~.—  a— - 
•  • 

5  sin  t  S  sin  t  _ 

/  . 


|  sin  t 


pb2Cg'dz 

7"  “  2*~ 

2m  sin  t  cos  t 


which  can  he  easily  integrated  to  give 


\  sin  t  ,  z0  +  I  sin  T  I  sin  t  w  z  +  S  sin  r 

- ; -  +  log  e  ;  .  .  6e  • 


zn  +  |  sin  t 


^  sin  t  z  +  |  sin  t 


|  sin  t 


pb 

.  2  2 
2m  sm  t  cos  t 


Cg  dz 


The  quantity  on  the  right-hand  side  of  equation  (12),  which  depends  only 
on  the  draft  (for  a  given  hull  shape,  weight,  and  trim),  will  henceforth 
he  abbreviated  as 


2  2  / 
2m  sin  t  cos  -r  2 


Cg  dz 


2U 
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M 


and  will  be  referred  to  as  the  generalized  draft.  An  alternate  form 
of  equation  (13),  in  terms  of  the  beam-loading  coefficient,  is 


k 


(14) 


where 


CA  -  (15) 

pb3 

In  the  remainder  of  this  paper  Cg  is  assumed  to  be  a  known  function 
of  draft  which  has  been  determined  either  experimentally  or  theoreti¬ 
cally.  Since  all  the  other  factors  in  equation  (l4)  are  known  con¬ 
stants,  except  the  quantity  m^/m.  which  was  used  for  converting  Cg 

to  Cg1  and  which  can  also  be  expressed  as  a  function  of  draft  as  indi¬ 
cated  in  appendix  A,  the  quantity  k  can  also  be  considered  as  a  known 
function  of  draft. 


In  order  to  compute  the  vertical  hydrodynamic  force  from  equa¬ 
tion  (8),  equation  (12)  must  be  expressed  in  terms  of  £  according  to 
equation  (9)  as  follows: 


+  logp 


I  tan  t 


-  1’ 


k£  tan  t 


i 


+  loge 


i 


|  tan  t 


^ |  tan  t 


(16) 


(The  quantity  1  has  been  arbitrarily  added  and  subtracted  from  the  left- 
hand  side  of  this  equation  for  later  convenience. )  In  order  to  facili¬ 
tate  computation,  a  function  \|r  is  defined  as  follows: 


\|r(co)  =  i  +  loge  CD-I 


(17) 


Equation  (l6)  can  be  expressed  in  terms  of  the  ^-function  as  (compare 
eqs.  (16)  and  (17)) 


*  - 


♦t- 


i 


=  k 


(l8a) 


1 1  tan  t/  tan  r. 


10 


MCA  TN  2814 


or 


* U-  -  *4 


-  k 


tan  t  /  VI  tan  T  / 


or,  in  terms  of  the  inverse  function  of  \]/,  as 


5 


=  + 


|  tan  t 


♦ 


tan  r  > 


(18b) 


(18c) 


f 1 

J  to  \  y 

i_  = 

Y  .  /  K 
\|  tan  t/ 

lo 

t0 

|  tan  t 


(I8d) 


The  constant  l0/i  tan  t  will  he  henceforth  abbreviated  as  e.  Since 
i  (which  equals  l0)  is  equal  to  £ Q  cot(70  +  t)  (see  fig.  1(b)),  e 
can  also  be  expressed  in  the  form 


=  tan(7o  +  t) 
6  tan  t  tan  T 


(19) 


It  is  noted  that  this  constant  e  is  related  to  the  impact  approach 
parameter  k  of  reference  3  by  the  relation  €  =  1  +  where 


K 


sin  t 
sin  70 


cos (7  o 


+  t) 


The  values  of  the  ^-function  necessary  for  the  use  of  equation  (l8d) 
are  given  in  table  I  and  in  figures  2(a)  and  2(b).  The  procedure  for 
using  figures  2(a)  and  2(b)  is  illustrated  in  figure  2(c).  For  example, 

to  compute  the  quantity  i)f1[t(e)  -  k] ,  first  find  6  on  the  u>-axis. 

The  corresponding  ordinate  is  i|r(e).  Then  go  down  on  the  t-axis  a  dis¬ 
tance  k  to  obtain  i|r(e)  -  k.  The  quantity  t“15'(e)  -  k]  is  then  the 
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point  on  the  oo-axis  corresponding  to  i|r(e)  -  k  on  the  ty-axis.  (For 
more  accurate  results  a  similar  procedure  can  be  followed  with  the  aid 
of  table  I.)  It  is  noted  that  the  inverse  ^-function  is  a  double-valued 
function.  The  physical  significance  of  this  double-valued  nature  for 
the  present  problem  is  as  follows:  In  each  impact  each  value  of  draft 
(and  of  k)  is  reached  twice,  once  while  the  seaplane  is  going  down  into 
the  water  and  once  while  it  is  coming  out  of  the  water.  Values  of 

f1  i(e)  -  k]  greater  than  unity  correspond  to  the  descending  motion 

and  those  less  than  unity  correspond  to  the  ascending  motion. 

The  vertical  hydrodynamic  force  is  now  obtained  as  a  function  of 
the  generalized  draft  by  combining  equations  (8),  (l8d),  and  (19)  to 
yield 


F  = 


2'  2  2 

pb  |  tan  t  , 


0-2 
2  sin  t 


,-1 


V(e)  -  k 


^  2---2 


cos 


iio  <- t)  „ .  C-i 


O  2 

2  cos  T 


t(0  -  k 


(20) 


The  result  may  also  be  expressed  as  a  dimensionless  lift  coefficient 


CT  = 


cos 


1  2  2 
A  pv0  b 


0 o  +  T)  .  ,j,-i 

“CB  V 


2 

COS  T 


t(e)  -  k 


(21) 


or  as  a  load  factor 


b^  2--  2 


F  pb  V _  cos 
z  _  o 


xw  mg 


2mg  cos^t 


(?o-iCb,  r+-i[t(e).k] 


(22) 


The  vertical  velocity  can  be  related  to  the  generalized  draft  by 
combining  equations  (9),  (l8c),  and  (19)  to  give 


V  1|jHO  -  k]  -  ! 


e  -  1 


(23) 
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The  rebound  (and  initial)  vertical  velocity  is  obtained  by  setting 
the  generalized  draft  equal  to  zero: 


It  should  be  noted  that  the  rebound  velocity  depends  only  on  the  impact 
parameter  e  and  is  not  affected  by  the  shape  or  mass  of  the  seaplane. 

The  maximum  draft  is  obtained  by  setting  the  vertical  velocity 
equal  to  zero  in  equation  (23)  so  that  i  af(6)  -  w]  =  1  or 
\|r(e)  "  kmax  =  ^(l)>  and  since  i|r(l)  =  0  (see  table  I), 

w  *  +(«>  («) 

where  the  relation  between  J?max  and  zmax  is  given  by  equation  (l4). 

The  relation  between  the  time  and  the  draft  is  given  by 

t  =  f  —  dz 
Jq  dz 

that  is,  a  plot  is  made  of  the  quantity  l/z  against  z  and  is  graphi¬ 
cally  integrated  to  obtain  the  time. 

A  proposed  procedure  for  computing  impact  loads  and  motions 
according  to  the  preceding  equations  is  given  in  appendix  B. 

Substantiation  of  Theory 

The  results  of  the  preceding  derivation  can  be  compared  with  experi¬ 
mental  impact  loads  and  motions  in  order  to  test  the  validity  of  these 
results  and  thus  of  the  assumptions  upon  which  the  derivation  is  based. 
For  the  non- chine -immersed  case  the  theory  of  the  present  paper  is  sub¬ 
stantially  the  same  as  the  theory  of  references  1  to  6 ;  hence,  the  agree¬ 
ment  between  theory  and  experiment  demonstrated  in  those  papers  for  both 
the  planing  and  the  impact  cases  constitutes  a  verification  of  the  theo¬ 
retical  equations  of  this  paper.  However,  for  impacts  which  involve  more 
than  a  slight  degree  of  chine  immersion,  such  as  those  of  a  rectangular 
flat  plate,  no  substantiation  of  the  theory  has  been  previously  given. 
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To  supply  such  a  substantiation,  experimental  and  calculated  impact 
loads  and  motions  for  a  rectangular  flat  plate  are  compared  in  figure  3- 
Experimental  time  histories  of  draft,  vertical  velocity,  and  hydro- 
dynamic  load  are  shown  for  five  landings.  The  planing  data  needed  for 
these  calculations  were  obtained  from  reference  7  (see  appendix  C  for 
details  of  computations)  and  the  corresponding  impact  data  were  obtained 
from  references  8  and  9-  The  good  agreement  shown  in  -figure  3  for  the 
calculated  and  experimental  loads  and  motions  substantiates  the  general 
validity  of  the  impact  equations.  The  small  discrepancies  that  do  exist 
are  clearly  within  the  limitations  of  accuracy  of  the  experimental  impact 
measurements . 


DISCUSSION  OF  PLANING  DATA  AND  ANALYSES 


In  view  of  the  preceding  substantiation  of  the  impact  equations  of 
this  paper,  the  remaining  problem  is  to  provide  reliable  means  of  pre¬ 
dicting  the  planing  characteristics  of  seaplanes.  The  available  infor¬ 
mation  on  this  subject  can,  for  convenience,  be  grouped  into  the  fol¬ 
lowing  subdivisions:  (l)  the  rectangular  flat  plate,  (2)  prismatic 
bodies  with  angles  of  dead  rise  greater  than  10°,  (3)  prismatic  bodies 
with  angles  of  dead  rise  between  0°  and  10°,  and  (4)  prismatic  bodies 
of  arbitrary  constant  cross  section.  These  various  subdivisions  will 
now  be  considered  in  detail. 


The  Rectangular  Flat  Plate 

A  large  quantity  of  flat-plate  planing  data  are  available  in  refer¬ 
ences  7  and  10  to  13 .  However,  much  of  this  information  is  incomplete 
or  was  obtained  at  such  low  planing  speeds  that  the  buoyant  forces  (which 
are  considered  negligible  during  landings)  are  significant,  and  these 
data  cannot  always  be  safely  extrapolated  to  find  the  high-speed 
(buoyancy-free)  planing  characteristics  needed  for  impact  calculations. 
Also,  these  available  data  are  usually  limited  to  relatively  low  trims 
or  small  ratios  of  wetted  length  to  beam.  Consequently,  it  is  doubtful 
that  the  available  data  are  adequate  to  cover  all  practical  landing 
conditions . 

The  available  analyses  of  flat-plate  planing  data,  notably  refer¬ 
ences  14  to  l6,  are  all  primarily  empirical.  They  are  consequently  use¬ 
ful  for  interpolating  the  experimental  data,  but  they  cannot  necessarily 
be  extrapolated  accurately  beyond  the  range  of  the  experimental  data. 

One  difficulty  in  using  these  empirical  analyses  is  the  fact  that  they 
are  restricted  to  the  determination  of  the  relation  between  the  planing 
coefficient  C-g  and  the  ratio  of  wetted  length  to  beam,  whereas  for 

impact  calculations  the  relation  between  C-g  and  the  draft  must  be 
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known.  In  the  case  of  the  flat  plate  the  water  rises  up  in  front  of 
the  plate  (see  fig.  1(d))  so  that  the  wetted  length  Xe  materially 
exceeds  the  length  of  the  model  below  the  undisturbed  water  surface  A.^, 

where  A,j  is  related  to  the  draft  by  the  equation  =  - - . 

a  u  b  sin  t 

Planing  data  showing  the  relation  between  Xg  and  A,^  are  plotted  in 

figure  4 ,  together  with  a  simple  empirical  equation  fitting  the  data. 
This  relation  appears  to  be  substantially  independent  of  the  trim  for 
the  range  of  data  available,  namely,  for  trims  below  18°.  For  higher 
trims,  according  to  impact  data  in  reference  9,  "the  wave  rise  A,e  -  A,^ 

increases  with  increase  of  trim.  (The  quantity  A.e  used  in  this  paper 
roughly  corresponds  to  the  quantity  A-p  in  reference  9-  ) 


Prismatic  Bodies  With  Angles  of  Dead  Rise  Greater  Than  10° 

The  problem  of  the  determination  of  the  planing  characteristics  of 
a  V-bottom  hull  where  the  chines  do  not  penetrate  the  water  surface  has 
been  treated  rather  extensively,  both  theoretically  and  experimentally, 
in  references  1,  3,  5,  and  6  for  angles  of  dead  rise  between  10°  and  40°. 
However,  when  the  chines  do  penetrate  the  water  surface  the  analysis 
becomes  more  difficult.  The  first  important  problem  is  to  determine 
when  chine  immersion  occurs  or,  more  specifically,  the  draft  at  which 
chine  immersion  occurs  zc^.  Figure  l(c)  shows  that  this  occurs  when 

the  splashed-up  water  crosses  the  chines,  or  the  draft  is 


H  cos  t  b  tan  3  cos  t 
Zch  R  ~  2lT 


(27) 


(in  fig.  1(c)  water  pile-up  in  front  of  the  hull  has  been  neglected,  or 
A^  =  A,e.)  The  wave-rise  factor  R  has  been  evaluated  in  different  ways 

by  several  investigators.  According  to  Wagner’s  two-dimensional  deri¬ 
vation  for  small  angles  of  dead  rise  (ref.  17),  R  =  ^.  This  value  of  R 

is  supported  by  the  theoretical  solutions  of  reference  18,  which  also  show 
that  R  decreases  with  increase  of  dead-rise  angle,  and  by  the  experi¬ 
mental  data  of  reference  15  for  angles  of  dead  rise  between  10°  and  30°. 
Some  full-scale  landing  data  from  a  flying  boat  with  20°  angle  of  dead 
rise  (ref.  19)  indicate  a  value  of  R  =  1.4.  The  impact  theory  of  refer¬ 
ence  3  uses  a  value  of  R  =  Qj-  -  l^tan  p  and  reference  6  uses  a  value 

of  R  =  1.  It  should  be  noted  that  the  last  two  values  of  R  do  not 
directly  correspond  to  the  physical  water  splash-up  but  are,  rather. 
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effective  values  of  splash-up.  The  first  value  of  R  was  proposed  from 
considerations  of  the  over-all  load  problem  and  the  second  value  resulted 
from  the  semiempirical  analysis  of  planing  data  given  in  reference  6.  In 
view  of  these  observations  it  appears  that  R  is  approximately  equal 
to  tt/2;  however,  further  research  on  this  subject  would  be  of  value. 

Various  analyses  have  been  made  for  the  planing  coefficients  of  the 
V-bottom  hull  with  chines  substantially  immersed,  notably  those  of  refer¬ 
ences  6,  15,  and  20.  Reference  6  presents  a  semiempirical  treatment  of 
the  problem  which  also  considers  the  effects  of  buoyancy  on  the  planing 
loads.  In  order  to  bring  this  treatment  into  agreement  with  the  extremely 
limited  data  available  at  that  time  (refs.  6,  11,  and  13)  it  was  necessary 
to  assume  that  the  wave-rise  factor  R  was  equal  to  1.  Recently,  more 
V-bottom  planing  data  have  been  published  (ref.  15),  and  from  a  consider¬ 
ation  of  these  data  as  well  as  from  the  consideration  of  some  impact  data 
(to  be  discussed  subsequently)  for  a  float  with  30°  angle  of  dead  rise, 
the  analysis  of  reference  6  still  appears  reasonable  for  its  two  end¬ 
point  conditions,  the  non-chine- immersed  case  and  the  deeply  chine- 
immersed  case,  but  it  is  somewhat  conservative  in  the  intermediate  range 
as  a  result  of  the  assumption  that  the  wave-rise  factor  is  unity.  In 
reference  15,  where  the  more  recent  data  are  used,  an  empirical  analysis 
restricted  to  the  case  of  the  chine-immersed  hull  is  given  which  corre¬ 
sponds  to  the  more  realistic  splash-up  factor  of  jt/2  and  which  is  in 
fair  agreement  with  the  available  experimental  data.  (A  somewhat  parallel 
but  less  detailed  treatment  is  given  in  reference  20. )  However,  since 
the  available  data  are  still  not  very  extensive  and  do  not  adequately 
cover  the  practical  limits  of  large  ratios  of  wetted  length  to  beam,  it 
is  not  yet  certain  that  the  empirical  formula  of  reference  15  can  be 
safely  extrapolated  to  large  ratios  of  wetted  length  to  beam.  There¬ 
fore,  a  need  apparently  exists  for  further  study  of  the  planing  charac¬ 
teristics  of  V-bottom  hulls  with  angles  of  dead  rise  greater  than  10° 
for  the  case  of  considerable  chine  immersion. 

As  a  test  of  the  relative  merits  of  the  two  planing  equations  of 
references  6  and  15 >  calculations  based  on  these  two  formulas  were  made 
for  the  landing  loads  on  a  prismatic  V-bottom  hull  having  an  angle  of 
dead  rise  of  30°  and  are  compared  with  experimental  impact  data  in  fig¬ 
ure  5.  (Since  the  analysis  of  reference  15  is  restricted  to  the  chine- 
immersed  case,  the  planing  coefficients  prior  to  chine  immersion  were  in 
both  cases  determined  from  the  equation  of  reference  6  as  is  discussed 
in  appendix  D. )  At  the  low  flight-path  angle  (yQ  =  2.2°)  of  figure  5(a), 
which  corresponds  to  a  landing  in  which  the  chines  become  immersed,  but 
not  deeply,  the  analysis  of  reference  15  gives  much  better  agreement  with 
the  experimental  impact  data  both  in  the  shape  of  the  load  time  history 
and  in  the  magnitudes.  This  better  agreement  is  due  to  the  fact  that 
reference  15  makes  use  of  a  more  realistic  estimate  of  the  draft  at  which 
chine  immersion  begins  than  does  reference  6.  From  an  examination  of  fig¬ 
ure  5(h),  which  gives  computed  and  experimental  data  for  a  large  flight- 
path  angle  (7 0  =  5-9°)  which  corresponds  to  substantial  chine  immersion. 
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both  analyses  give  reasonable  qualitative  representations  of  the  load 
time  history,  reference  6  being  somewhat  conservative  and  reference  15 
being  somewhat  unconservative.  Additional  calculated  and  experimental 
impact  load  data  shown  in  figure  5(c)  substantiate  these  observations. 
In  figure  5(d)  are  shown  calculated  and  experimental  maximum  drafts. 
These  indicate  the  same  conclusions  as  the  load  comparison;  the  maximum 
draft  is  underestimated  by  reference  6  (which  is  the  same  as  over¬ 
estimating  the  load)  and  reference  15  tends  to  overestimate  the  draft. 

To  summarize  the  preceding  discussion  of  the  planing  coefficients 
for  angles  of  dead  rise  greater  than  10°,  the  use  of  the  equations  of 
references  6  and  15  apparently  results  in  a  reasonable  calculation  of 
the  landing  loads,  reference  6  giving  somewhat  conservative  results  and 
reference  15  giving  somewhat  unconservative  results.  However,  for 
greater  accuracy  further  planing  research  seems  to  be  needed. 


Prismatic  Bodies  With  Angles  of  Dead  Rise  Between  0°  and  10° 

As  yet  no  experimental  data  are  available  to  test  the  existing  anal¬ 
yses  for  angles  of  dead  rise  between  0°  and  10  .  The  analysis  of  refer¬ 
ence  6  is  known  to  be  inadequate  for  the  case  of  0°  angle  of  dead  rise; 
therefore,  it  is  as  yet  somewhat  doubtful  that  this  equation  can  be  used 
in  this  range  (0°  <  p  <  10°).  The  equation  of  reference  15,  by  the 
nature  of  its  derivation,  can  be  reasonably  assumed  to  apply  for  dead- 
rise  angles  in  this  range.  However,  one  difficulty  arises  in  the  use  of 
this  equation  in  this  range.  The  results  of  reference  15  are  essentially 
expressed  in  such  a  form  that  the  planing  coefficient  is  related  to  the 
wetted  length  at  the  keel  (Xe  in  fig.  1(d))  rather  than  to  the  draft. 

For  angles  of  dead  rise  greater  than  10°,  reference  15  states  that  the 
relation  between  the  draft  and  the  wetted  length  at  the  keel  is  approxi¬ 
mately  equal  to  the  sine  of  the  trim  ^Xe  =  However,  for 

the  case  of  the  flat  plate  (0°  angle  of  dead  rise)  this  relation  was 
previously  shown  to  be  considerably  more  complicated  because  of  the  wave 
rise  in  front  of  the  plate  (see  fig.  1(d)).  Thus,  if  reliable  computa¬ 
tions  are  to  be  made  for  the  landing  loads  for  V-bottom  surfaces  with 
angles  of  dead  rise  between  0°  and  10°,  a  need  exists  for  an  investiga¬ 
tion  of  the  wave-rise  effects  in  this  range. 


Prismatic  Bodies  of  Arbitrary  Constant  Cross  Section 

For  prismatic  hulls  of  arbitrary  constant  cross  section,  little 
experimental  information  is  available  on  planing  characteristics  (ref.  11). 
However,  reference  2  presents  a  theoretical  means  of  analysis  which  is 
applicable  to  hulls  of  constant  cross  section  without  reflex  chines  or 
chine  strips,  and  which  has  been  partly  substantiated  by  some  experimental 
impact  data  in  the  same  reference. 
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CONCLUDING  REMARKS 


The  equations  for  the  landing  impact  of  a  seaplane  with  a  hull  of 
arbitrary  constant  cross  section  have  been  presented  in  such  a  form 
that  the  landing  loads  and  motions  are  expressed  in  terms  of  the  planing 
properties  of  the  seaplane.  Comparisons  of  calculated  and  experimental 
loads  and  motions  for  the  landing  of  a  rectangular  flat  plate  have  shown 
that  this  analysis  is  valid  for  a  greater  range  of  conditions  than  was 
covered  by  previous  papers.  A  survey  and  evaluation  of  the  available 
information  on  seaplane  planing  characteristics  which  is  suitable  for 
use  with  the  analysis  of  this  paper  has  indicated  a  need  for  further 
research  on  these  characteristics. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va. ,  July  7*  1952 
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APPENDIX  A 

VIRTUAL  MASS 


In  reference  2  an  approximate  expression  is  given  for  the  virtual 
mass  of  a  body  of  arbitrary  constant  cross  section  prior  to  chine  immer¬ 
sion.  This  equation,  converted  to  the  terminology  of  the  present  paper, 
is 


m 


2C, 


sm  T  COS  T 


'B 


d£ 


(Al) 


This  relation,  substituted  into  equation  (l4),  gives 


k  = 


2  2 

2C^  sin  t  cos  t 


CB 


1  + 


s/b 


.  P  P 

2C^  sin^T  cos  t 


CB  *§ 


loge 


2C, 


1 


2 

sin  t  cos 


2 

T 


(A2) 


prior  to  chine  immersion. 

Subsequent  to  chine  immersion  the  virtual  mass  m,r  associated 

"ch 

with  the  wetted  area  forward  of  the  intersection  of  the  chines  and  the 
water  surface  can  be  computed  from  equation  (Al).  The  virtual  mass  of 
the  region  aft  of  this  intersection  can  be  approximated  by  a  semicylinder 
of  water  circumscribed  about  the  beam  and  corrected  for  the  finite  aspect 
ratio.  The  resulting  equation  for  the  virtual  mass  is 


pjrb“ 


n»"^ch+  8 


XCh9(A) 


(z  ^  z 


ch 


)  (A3) 
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or 


(z  2-  zch)  (A4) 


where  cp(A)  is  the  aspect-ratio  correction.  For  cases  where  the  water 
does  not  splash  up  in  front  of  the  model,  the  wetted  chine  length  X-ch 

(see  fig.  1(d))  is  equal  to  - —  and  equation  (A4)  becomes 

b  sin  t 


% 

m 


on  .  2  2 

2C^  sin  t  cos  t 


rcp(A)  /z  _  fch] 
8Ca  sin  b  / 


(A5) 


For  the  flat  plate  the  wetted  chine  length  is  the  total  wetted  length 
(lch  =  ^e)  and  equation  (A4)  becomes 

T  =  7=-  MM  (AS) 

COA 

For  arbitrary  shapes  with  splash-up  in  front,  the  use  of  equation  (A5) 
is  suggested.  The  error  introduced  by  neglecting  the  splash-up  for  this 
purpose  will  usually  be  small  or  negligible. 

According  to  Pabst  (refs.  21  and  22)  the  virtual  mass  of  a  rectan¬ 
gular  flat  plate  of  aspect  ratio  A,  as  determined  by  vibration  tests, 
is  equal  to  that  of  a  circumscribed  cylinder  times  the  quantity  (p(A), 
where 


9(A) 


(0  <  A  <  oo ) 


(AT) 


which  can  be  approximated  as 


qp(A)  =  i 


_i_ 

2A 


(A  >  1.5) 
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These  relations  are  suggested  for  use  in  the  computations  of  the  aspect- 
ratio  correction  for  models  of  arbitrary  shape,  provided  that  the  aspect 
ratio  of  the  model  can  he  computed  as  a  function  of  draft  without  too 
much  difficulty.  If  this  procedure  is  difficult,  however,  the  aspect- 
ratio  correction  may  be  considered  approximately  equal  to  unity: 

9(A)  ~  1  (A9) 

For  cases  where  A  is  reasonably  large,  no  significant  error  should 
result  from  the  substitution  of  equation  (A9)  into  equation  (A5). 

The  hydrodynamic  aspect  ratio  of  a  model  is,  by  definition,  equal 
to  the  quotient  of  the  square  of  the  keel  wetted  length  and  the  pro¬ 
jected  wetted  area  normal  to  the  keel.  Thus,  for  a  rectangular  flat 
plate. 


A  =  Xe  (A10) 

and,  for  a  wedge  of  dead-rise  angle  f3  with  no  forward  splash-up  of 
water, 

(tf  CiiC  T 

tan  (3  cos  t 

4r 


z 

b 


(All) 
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APPENDIX  B 

PROPOSED  PROCEDURE  FOR  COMPUTATIONS  OF  IMPACT  LOADS  AND  MOTIONS 


In  order  to  utilize  the  equations  of  this  paper  to  compute  the 
water  loads  and  motions  of  a  seaplane,  the  following  procedure  is 
suggested: 

(1)  Select  values  of  C^,  t,  yQf  and  VQ. 

(2)  Make  a  plot  of  CB  against  z/b.  This  plot  may  be  obtained 

from  experimental  planing  data  or  by  utilizing  theoretical  planing 
e  quat ions . 

(3)  Graphically  or  analytically  integrate  the  CB  curve  to  obtain 

rz/b 

the  quantity  J  CB  d|  as  a  function  of  draft. 


(4)  Compute  mw/m  as  a  function  of  draft  according  to  equations  (Al) 
and  (A5).  (For  slightly  better  accuracy  in  the  case  of  the  rectangular 
flat  plate,  use  equation  (A6),  together  with  the  empirical  curve  in  fig¬ 
ure  4,  instead. ) 

(5)  Compute  CB’  as  a  function  of  draft  according  to  equation  (7). 

(6)  Compute  k  as  a  function  of  draft  according  to  equation  (l4). 

(7)  Compute  e  from  equation  (19)  and  then  \|r(e)  from  table  I 
or  figure  2. 

(8)  Select  a  series  of  values  of  k  between  0  and  t(e). 

(9)  For  each  value  of  k  (the  generalized  draft),  compute 

i  1[t(6)  _  kj  as  indicated  in  figure  2(c),  with  the  aid  of  table  I  or 

figure  2(a)  or  2(b).  ^The  larger  value  of  3  obtained 

corresponds  to  the  seaplane  descending  into  the  water  and  the  smaller 
value  corresponds  to  the  seaplane  rising  out  of  the  water. j 

(10)  Compute  the  vertical  force  from  equation  (20),  (21),  or  (22). 
Plot  this  quantity  against  draft.  The  maximum  load  can  then  be  obtained 
from  this  curve. 


22  NACA  TN  2814 

(11)  Obtain  the  vertical  velocity  as  a  function  of  draft  from 
equation  (23). 

(12)  Obtain  the  relation  between  draft  and  time  by  graphical  inte¬ 
gration  of  equation  (26). 
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APPENDIX  C 

DETAILS  OF  COMPUTATIONS  FOR  FLAT  PLATE 


The  impact  loads  and  motions  shown  in  figure  3  for  the  rectangular 
flat  plate  were  computed  according  to  the  procedure  of  appendix  B.  The 
relations  between  Cg  and  z/b  needed  for  step  (2)  of  that  procedure 

were  obtained  as  follows.  The  relation  between  the  wetted  length  Xe 
and  the  planing  coefficient  Cg  was  obtained  from  a  cross  plot  of 

Locke's  high-speed  planing  data  from  reference  7  (see  fig.  6).  The 
relation  between  the  draft  and  the  wetted  length  was  obtained  from  fig¬ 
ure  4,  which  is  a  plot  of  data  obtained  by  Sottorf  (ref.  10)  and 
Sambraus  (ref.  12).  The  empirical  formula  shown  in  this  figure  was  used 
for  all  computations.  The  relation  between  C-g  and  z/b  was  obtained 

by  combining  figure  6  with  the  empirical  formula  from  figure  4. 
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APPENDIX  D 

DETAILS  OF  COMPUTATIONS  FOR  30°  ANGLE  OF  DEAD  RISE 


The  impact  loads  and  motions  shown  in  figure  5  for  a  model  with 
30°  angle  of  dead  rise  were  computed  according  to  the  procedure  of 
appendix  B.  The  relation  between  CB  and  z/b  needed  for  step  (2)  of 

that  procedure  was  obtained  from  the  following  equations  taken  from 
references  6  and  15: 

According  to  reference  6  (if  buoyancy  be  neglected). 


CB  =  1.42*  cot2p  sin  r(l  -  (z  <  zch) 


CB  =  0. 355*  sin  r  cos^t 


/z  tan  p\ 

Sin  T - r— - 

, _ \b  4  I 

-4  — 


Hi 


>  (Dl) 


„  .  Jz  tan  p  cos  t\ 

C8  sin  t  cos  t^— - - - J 


(z  >  zch) 


R  =  1 


where  Ca  is  a  function  of  dead-rise  angle  tabulated  in  reference  6 
(0.750  for  3  =  30°). 

According  to  reference  15,  for  the  chine -immersed  case  only  (if 
buoyancy  be  neglected), 


CB  =  C-g^  —  00653Cg 


0.6 


r  (z  >  zch) 


R  =  5 


(D2) 
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where 


cBf  =  O.OlPOh/^T1-1 

^-m  =  ^(^e  +  ^ch) 

>  -  z 

*“  -u 

b  sin  t 


Xch  -  xe 


tan  3 
ir  tan  t 


Since  equation  (D2)  is  not  valid  for  the  case  of  no  chine  immer¬ 
sion,  in  order  to  use  this  equation  for  impact  computations  equation  (Dl) 
was  used  to  predict  the  planing  coefficient  prior  to  chine  immersion. 
Chine  immersion  was  taken  to  occur  at  the  draft  at  which  equations  (Dl) 
and  (D2)  predicted  the  same  planing  coefficient.  This  procedure  for 
determining  chine  immersion  was  chosen  only  because  for  the  specific 
case  considered  (p  =  30°,  t  =  15°)  this  value  of  draft  corresponded  to 
a  wave-rise  factor  of  roughly  it/2,  which  is  in  agreement  with  the  value 
predicted  by  equation  (D2). 
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Figure  1.-  Geometrical  relations 
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Figure  2.-  Variation  of  the  -function.  i (cjo)  =  —  +  log  cd  -  1. 
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Figure  3.-  Comparison  of  calculated  and  experimental  loads  and  motions 
for  a  rectangular  flat  plate  with  a  beam-loading  coefficient  of  l8.8 
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O  Experimental  data  obtained  from  tests 

_ Impact  computations  based  on  analysis 

_ Impact  computations  based  on  analysis 


described  in  reference 
of  reference  6 
of  reference  15 


23 


(c)  Maximum  loads. 


(d)  Maximum  drafts. 


Figure  5.-  Comparison  of  calculated  and  experimental  loads  and  drafts 
Sfor  a  prismatic  model  with  an  angle  of  dead  rise  of  30  ,  a  beam- 
loading  coefficient  of  l8.8,  and  a  trim  of  15  • 
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